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ABSTRACT 

Models of disc galaxies which invoke viscosity-driven radial flows have long been 
known to provide a natural explanation for the origin of stellar exponential discs, 
under the assumption that the star formation and viscous timescales are comparable. 
We present models which invoke simultaneous star formation, viscous redistribution 
of gas and cosmologically-motivated gaseous infall and explore the predictions such 
models make for the scale length evolution and radial star formation history of galactic 
stellar discs. While the inclusion of viscous flows is essential for ensuring that the 
stellar disc is always exponential over a significant range in radius, we find that such 
flows play essentially no role in determining the evolution of the disc scale length. In 
models in which the main infall phase precedes the onset of star formation and viscous 
evolution, we find the exponential scale length to be rather invariant with time, with 
the bulk of the disc stars at all radii out to ~ 5 scale lengths being relatively old (ie. 
ages <; 6 — 8 Gyr for an assumed disc age of 11 Gyr). On the other hand, models 
in which star formation/ viscous evolution and infall occur concurrently result in a 
smoothly increasing scale length with time, reflecting the mean angular momentum 
of material which has fallen in at any given epoch. The disc stellar populations in 
these models are predominantly young (ie. ages <; 5 Gyr) beyond a few scale lengths. 
In both cases, viscous flows are entirely responsible for transporting material to very 
large radii. Our results are robust for a range of currently popular star formation laws 
and infall prescriptions. We discuss existing observational constraints on these models 
from studies of both local and moderate redshift disc galaxies. In particular, a good 
agreement is found between the solar neighbourhood star formation history predicted 
by our infall model and the recent observational determination of this quantity by 
Rocha-Pinto et al (2000). 

Key words: galaxies: formation - galaxies: evolution - galaxies: structure - galaxies: 
spiral - galaxies: fundamental parameters 



1 INTRODUCTION 

One of the most notable properties of disc galaxies is that the 
surface brightness profiles of their stellar discs are remark- 
ably close to exponential distributions over a large dynamic 
range in radii (eg. de Vaucouleurs 1959; Freeman 1970; de 
Jong 1996a). This has traditionally been held to imply an 
exponential surface density profile for gas infalling to the 
galactic plane, which in turn places particular constraints on 
the detailed angular momentum distribution of baryonic ma- 
terial in the protogalaxy. Gunn (1982), for example, showed 
that, under the assumption of angular momentum conserva- 
tion for each gas element under infall, a uniformly rotating 
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uniform sphere collapses to an approximately exponential 
profile, however the resulting function is only exponential 
over ~ 3 scale lengths compared to the more impressive 4-6 
scale lengths exhibited by real galaxies (eg. de Jong 1996a). 
Calculations with other initial halo density profiles and an- 
gular momentum distributions produce similar results (eg. 
van der Kruit 1987; Dalcanton, Spergel & Summers 1997; 
Bullock et al 2000). 

The exponential profile is so ubiquitous, and its appear- 
ance when plotted as magnitude versus radius so deceptively 
simple, that one may be tempted to regard it as a natural 
or obvious outcome of gravitational collapse. In fact, gravi- 
tational collapse calculations in other astronomical contexts 
never give rise to exponential disc profiles (Cassen & Moos- 
man 1981; Terebey, Shu & Cassen 1984), with the more usual 
outcome being profiles of approximately power law form. In 
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the galactic context, the difficuhy of producing exponential 
profiles over a large number of scale lengths through angular 
momentum conserving collapse from reasonable initial con- 
ditions has recently been highlighted by Efstathiou (2000): 
the outer regions of galactic discs contain so little mass, that 
for power law halo density profiles, the outer disc material 
must originate from a small range of radii in the halo. On 
the other hand, these regions contain a large dynamic range 
of angular momentum, so that these conditions can only be 
met simultaneously if the angular nionicntum profile at the 
halo suffers an abrupt, and physically implausible, upturn 
at some radius (see Figure 6(b) of Efstathiou 2000). 

Exponential profiles have the unique feature that the 
scale length on which quantities change is the same at all 
radii. In the case of a contrifugally-supported disc, the ex- 
ponential scale length is determined (in a fixed potential) 
purely by the average angular momentum per unit mass 
for the entire disc. Such a system thus has the remarkable 
property that material over a significant range in radii (and 
with a correspondingly significant range in specific angular 
momenta) is somehow 'informed' about the average specific 
angular momentum of the disc as a whole. The key question 
is whether this property is a result of the formation pro- 
cess of the disc, or somehow acquired through subsequent 
evolution. 

Disc formation models have moved on from the picture 
of smooth dissipational collapse within a dark halo poten- 
tial, and recent numerical calculations have focused on the 
more complex formation histories resulting from hierarchical 
growth (eg. KaufFmann, White & Guiderdoni 1993; Cole et 
al 1994; Navarro, Frenk & White 1995). A notable outcome 
of such calculations is the strong redistribution of angular 
momentum during collapse. Attention has so far mainly fo- 
cused on the net transfer of angular momentum from bary- 
onic to dark matter and the consequent implications for 
the overall sizes of galactic discs (eg. Navarro & Bonz 1991; 
Navarro & Steinmetz 1997; Weil, Eke & Efstathiou 1998). 
A related issue is, however, the redistribution of angular 
momentum within the baryonic material and its effect on 
the resulting disc profile. The large dynamic range in size 
scales that have to be captured in such simulations, and the 
difficulties in controlling the effective viscosity in SPH cal- 
culations that are close to the limit of adequate resolution, 
mean that it is not yet clear whether discs that are exponen- 
tial over a number of scale lengths are indeed the expected 
outcome. 

Observational evidence against a picture of exponential 
infall followed by in situ star formation is provided by the 
radial profiles of gas in disc galaxies. The gas surface density 
profiles are generally observed to be much flatter then the 
stellar profiles. In the inner regions of disc galaxies, the gas is 
the minority component and represents the residue following 
star formation: it is not difficult, by appropriate adjustment 
of the star formation law, to generate gas profiles within 
the optical disc that are compatible with observations using 
an initially exponential gas profile (eg. Boissier & Prantzos 
1999). In the outer regions, however, the gas is the dominant 
component and its profile should presumably reflect its infall 
distribution. It seems highly suspicious, in the context of in 
situ models, that the gas profiles at large radii are so different 
from the exponential profiles that characterise the stars at 
smaller radii: if the stellar exponential profile derives from 



exponential infall of gas, surely one would expect to see that 
profile continued in the gas not yet incorporated in staxs? 

All of these difficulties are circumvented by the sim- 
ple model first proposed by Lin and Pringle (1987a). These 
authors envisaged that gas in the disc not only turns into 
stars (on a timescale t*), but is also subject to some form of 
effective viscosity that gives rise to angular momentum re- 
distribution and radial flows. The dominant physical process 
responsible for the viscosity is poorly constrained, but could 
plausibly result from cloud-cloud collisions and/or gravita- 
tional instabilities in the disc. If viscous redistribution oc- 
curs on a timescale tv, then the resultant stellar profile is 
exponential for a wide range of initially less centrally con- 
centrated gas profiles, provided only that ti, ^ t,. This result 
can be understood, in the context of the above discussion, 
inasmuch as viscous processes in the disc (provided t* is not 
much less than t^) allow material at different specific angu- 
lar momenta to communicate with each other, and thus to 
'inform' themselves of the average specific angular momen- 
tum of disc material, which in turn sets the disc scale length. 
On the other hand, t» cannot be much greater than t„ be- 
cause in that case, viscous evolution would run away too far 
towards its ultimate ondpoint (i.e. 'all the mass at the ori- 
gin, all the angular momentum at infinity', Lynden-Bell and 
Pringle (1974)) to be compatible with the observed profiles 
of disc galaxies. Fortunately it turns out that the generation 
of exponential profiles over a number of scale lengths does 
not require particularly fine tuning between t» and t^,, and 
that satisfactory fits are obtained even if the ratio t,/ti, de- 
viates from unity by half an order of magnitude, or is a mild 
function of radius (Clarke 1989). 

The viscous model thus relieves cosmological models 
of the burden of generating the 'right' angular momentum 
profile (though, by the same token, it also erases informar 
tion about the profile on infall, and thus weakens the role 
of gaseous and stellar disc profiles as cosmological probes). 
It also predicts that the radial gas surface density distri- 
bution beyond the optical radius should be rather flat, in 
accordance with observations. The viscous model was able 
to satisfy the various observational constraints available a 
decade ago (mainly surface density profllos and chemical 
constraints imposed by the Milky Way), although it was 
not manifestly superior to conventional in situ models in 
this respect (Clarke 1989, Yoshii & Sommer-Larsen 1989). 
Then, as now, the main drawback to viscous models is the 
lack of a detailed theory linking to t,^, although it is not 
difficult to invoke physical processes (e.g. self-gravitational 
instability of the disc) that one might expect to simultane- 
ously promote star formation and angular momentum redis- 
tribution in the disc (Lin & Pringle 1987b; Olivier, Primack 
& Blumcnthal 1991). 

A wealth of new observational data pertaining to the 
evolution of disc galaxies has emerged in recent years from 
both high redshift galaxy surveys - which permit a direct 
study of how galactic properties evolve with epoch - as well 
as detailed studies of the fossil stellar populations in our 
own Milky Way and other local systems. Given that the 
general arguments in favour of viscous models appear at 
least as strong now as they did a decade ago, it would seem 
timely to revisit this model and investigate its predictions 
in the light of new and forthcoming observations. In this 
paper, we explore the structural evolution and radial growth 
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of stellar discs subject to simultaneous star formation and 
viscous evolution by considering models in which the initial 
gas disc is pre-assembled as well as models in which the 
gas disc is assembled over an extended period due to the 
continued infall of material from the protogalactic halo. In 
both cases, the gaseous and stellar discs are naturally built 
up from the inside-out, though the mechanism is different 
in each case: in pre-assembled models, gas is transferred to 
larger radii only as a result of viscous torques in the disc, 
whereas in the case of continued infall, the dominant cause 
of the inside-out growth is the increasing specific angular 
momentum of infalling material at later times (whether or 
not viscous flows are also important). 

We do not address in this paper how well viscous models 
can be made to reproduce in detail the present-day proper- 
ties of the Milky Way and other nearby spirals; for a recent 
discussion of this aspect see Thon & Meusinger (1998). We 
also bypass the issue of how one sets up the conditions in 
the collapsing gas that will reproduce the correct absolute 
values of the disc scale lengths. The resultant scale length in 
viscous evolution models is set purely by the average specific 
angular specific angular momentum of the gas infallen at a 
given time hence the problem of forming disc systems with 
large scale lengths (as discussed by Mo et al 1998) is one 
of cosmological initial conditions and not one that viscous 
evolution can alleviate or affect in any way. In the present 
paper, we simply assume that the forming disc has access 
to material of the requisite high angular momentum. We 
describe in Section ^ the model and free parameters and 
present our results in Section ^. Section ^ discusses model 
constraints from existing and future observational data on 
both local and distant spirals. 



2 MODEL PARAMETERS 

We envisage a general scenario for disc galaxy formation in 
which infalling material gradually builds up the disc over a 
time interval t = to t = tinfaii and (simultaneous) star 
formation and viscous redistribution occur over the inter- 
val t = tsf (< tinfall) to t ^ t„ow (> tinfall)- The CaSC 

tsf = tinfaii corresponds to the case of star formation in a 
pre-assembled gas disc whereas t now ^ iinfall ^sf COrrC- 
sponds to the extreme infall case where gas is accreted up 
until the present epoch. 

2.1 Infall Calculation 

We calculate the surface density profile with which material 
arrives in the disc plane by assuming an idealised (spheri- 
cally symmetric) density distribution and angular velocity 
distribution for baryonic material in the protogalaxy. We 
then map this material down into the disc plane by invoking 
detailed angular momentum conservation, neglecting both 
the self-gravity of the baryonic material and the response 
of the halo to the infalling baryons. We justify this rigid 
halo approximation on the grounds that the initial distri- 
bution of baryons in the halo is poorly known (and likely 
to be considerably more complex than the smooth spheri- 
cally symmetric distribution invoked here) and also because 
in any case, the behaviour of the viscous evolution is only 
weakly dependent on the shape of the initial distribution (as 



opposed to its total mass and angular momentum). We cal- 
culate the gas infall rate (for t < tinfaii) by requiring that 
all material infalling from within a spherical radius ri„f{t) 
has reached the disc plane by a time t. This infall radius is 
fixed by the requirement that the free fall time within ri„/ is 
much less than the instantaneous cosmic shear time {H~^) 
- or, equivalently, that the mean enclosed density is much 
larger (« 200) than the mean density of the Universe at 
that epoch. As successive shells of matter are 'released' and 
allowed to fall in, we map them onto the disc plane (with- 
out applying the (small) correction for the finite infall time 
of each shell). We do not explicitly consider the gas cooling 
time in our model; numerical simulations suggest that the 
cooling radius tracks the virial radius for Milky Way-type 
galaxies hence that all gas within ri„f is indeed able to infall 
on a free fall time (Somerville & Primack 1999). Our default 
halo model assumes that the baryons initially follow an 
density distribution (ie. isothermal sphere), so that for an 
Einstein-de Sitter cosmology, the total infall rate onto the 
galaxy is constant in time for t < tinfaii- At any given ra- 
dius however, the infall history is characterised by an initial 
lag of variable length (corresponding to the time required 
for a shell of high enough angular momentum to fall in), a 
peak and then a roughly exponential decline. This prescrip- 
tion thus produces a crude on/off switch on the infall, so 
that, if tinfaii is comparable to tnom, it implies considerably 
more infall at recent times and at large radii than is usually 
assumed in galactic infall models, where the infall is arbi- 
trarily assumed to decline exponentially in both time and 
radius (eg. Boissier & Prantzos 1999; Chiappini et al 1997) 
in order to mimic inside-out growth. We motivate the pre- 
scription employed here by the wish to explore the response 
of viscous discs to extreme infall conditions. 

Finally, in order to prescribe the surface density distri- 
bution generated by each shell as it falls in, we must assign 
an initial angular momentum distribution for the baryons 
in the halo. The angular momentum of the dark matter and 
baryons is presumably acquired through tidal torques from 
neighbouring structures (Peebles 1969) and thus the average 
angular momentum of infalling material (whether the infall 
be smooth or discrete) is expected to increase with time. We 
assume the halo baryons initially rotate rigidly on spheres 
with angular velocity oc r~^ , which implies that the mean 
value of the spin parameter A for each shell is constant. For 
a power law angular velocity law, the surface density distri- 
bution generated by each shell is self-similar. We scale the 
overall dimensions and normalisation of the system to match 
the observed profile of the Milky Way. As stated before, our 
calculation thus sidesteps the difficulty of how, given the 
usual assumptions about the value of A and the size scale 
for collapse as a function of redshift, it is possible to produce 
large. Milky Way-scale, discs at redshifts greater than ~ 1 
(eg. Mo et al 1998). 

Figure ^ shows the final gas surface density profile at 
the completion of the infall phase (hence the initial starting 
point in the case of a pre-assembled gas disc). Beyond the 
central region, the gas surface density slowly declines and 
then turns over abruptly at the maximum infall radius, rinit- 
The value of rinit /^d depends primarily on the initial halo 
profile and rotation law, and to a lesser extent on the star 
formation and viscosity prescriptions. If the disc surface den- 
sity decreases outwards throughout the disc, rinit > 3 Rd', 
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Figure 1. The final gas surface density profile at the completion 
of the infall phase (solid line) - before star formation switches on 
— for a singular isothermal halo. The bottom radial coordinate is 
normalised to the scale length, Rgas derived from an exponential 
fit to the intermediate region of the profile (dashed line). Note 
that an exponential provides a good fit to the profile only over 
~ 2.5 scalelengths. The top radial coordinate is normalised to the 
scale length of the stellar disc which ultimately forms from this 
gas distribution, for the case of a star formation law oc JlgasR~^ ■ 
For the halo profile and rotation law adopted here, ri^it ~ 5 Rd- 



the absolute minimum value of rinu, which corresponds to 
an annulus of matter at this location, is 2 Rd- This lat- 
ter configuration does not viscously evolve to a smooth ex- 
ponential profile however. For the flat rotation curve and 
halo profile assumed here (and neglecting disc self-gravity), 
rinit ~ 5 Rd- An exponential fit (dashed line) provides a 
good description of the radial behaviour of the initial gas 
profile over ~ 2.5 scale lengths. 



2.2 Disc Evolution 

The evolution of gas surface density, Egas {r, i) in a disc that 
is subject to simultaneous viscous redistribution and star 
formation is described by: 



dt 



ld_ 

r dr 



d/dr[uEgasr^{dn/dr)] 



^gas 

u 



(1) 



where r, Q., v and t* represent the galactocentric radius, 
angular velocity, kinematic viscosity and local star forma- 
tion timescale. The evolution of the stellar surface density 
S*(r, t) is then given by: 



(2) 



The gas surface density in these models is assumed to 
be the sum of the neutral and molecular components. We 
calculate the evolution of the gaseous and stellar components 
of the disc, for times t > tsf, assuming that the gas evolves 
in the fixed potential of a halo with flat rotation law (i.e. 
SI oc r~^). This simpliflcation will not be correct for the 
very inner regions of the disc, and so we do not attempt 



any detailed comparison of our results with structures in 
the central parts of disc galaxies. 

In order to obtain a resulting exponential stellar pro- 
file, we follow Lin & Pringle (1987a) and fix the relationship 
between the star formation and viscosity prescriptions such 
that the characteristic timescale for material to move a frac- 
tional radius of order unity is comparable with the timescale 
on which it is converted into stars (i.e. t« ~ tv)- With this 
constraint, and given an assumed rotation law, the only free 
parameter is the star formation prescription. Our default 
star formation law relates the star formation rate per unit 
area with the gas column density, Sgas, as Y,gasR~^ , which 
corresponds to an n = 1 Schmidt law, modulated by the 
local angular frequency (eg. Wyse & Silk 1989). A prescrip- 
tion of this sort (with power-law n — 1.5) has recently been 
shown to be in good agreement with the observed radial star 
formation rate in the Milky Way (Boissier & Prantzos 1999). 
We also consider a model with a dependence of T.sfr oc Sga^s 
as well as a model in which the star formation law steepens 
consderably in the outer disc, Tjsfb.{R > 3Rd) oc T,gas/R^- 
While the former has been shown to provide a reasonable 
description of the star formation rates observed in the inner 
regions of galactic discs, the two component law may provide 
a better description of the steep declines in star formation 
rate observed at large radii, perhaps resulting from the flar- 
ing of the gas disc in these parts (eg. Ferguson et al 1998). 
We stress that in this paper we are seeking generic proper- 
ties of viscous models, and their inter-relation with contin- 
ued infall, that do not depend strongly on the star formation 
prescription employed. While varying the star formation pre- 
scription also requires changing the form of the viscosity law 
so that the approximate equality in the respective timescales 
is preserved, the dominant viscosity-producing mechanism 
in discs is so poorly constrained that this does not create 
any obvious physical problem^ 

We compute the evolution of a viscous star forming disc 
subject to infall using a standard first order explicit scheme 
with 200 gridpoints, equispaced in R with Rout/Rin ~ 150. 
Zero torque boundary conditions are applied at the inner 
and outer edge of the grid. In practice, we prescribe tsf and 
tinfaii and then pursue the calculation until a time t„ow at 
which the system is deemed to resemble the Milky Way - 
specifically, we require that the gas fraction at a radius of 
2.7 exponential scale lengths^ is ~ 0.4. 



3 RESULTS 

3.1 The Case of a Pre- Assembled Gas Disk 

We first consider the results of models in which star forma- 
tion (hence viscous redistribution) switches on after the ini- 
tial gas disc is already in place - ie. pre-assembled. This cor- 
responds to the scenario where all or most of the infall phase 
is complete before the onset of star formation {tsf > tinfaii)- 
Figure ^ shows snapshots of the gaseous and stellar surface 

t For a given star formation law of the form Yjsfr oc Sg^j^-R"'', 

the corresponding viscosity prescription is oc Ti'^as R^~^ 

t We adopt Rq = 8.0 kpc (Reid 1993) and Ri3=3.0 kpc (Sackett 

1997). 
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Figure 2. Snapshots of the gaseous (left) and stellar (right) surface density distributions at various fractional times (1/3, 2/3 and 1) of 
the star formation duration (which terminates at tnom) for the case of a pre-assembled gas disc. The radial coordinate is normalised to 
the stellar exponential scale length measured at the end of the model. The final solutions are those with the highest stellar profile at all 
radii and the lowest inner gas profile. The initial gas disc in these models extends to approximately 5 (final) exponential scale lengths. 
Different star formation laws are shown in top, middle and bottom panels. 



density distributions at various fractional times since the on- 
set of star formation for three different star formation pre- 
scriptions (EsFR oc T,gas/R, SsFi? oc Eg^f^ and a 'two com- 
ponent' law described by T.sfr oc Tigas/R", where n = 1 
out to 3 scale lengths and n = 2 beyond that). The pro- 
files are normalised by matching the final (ie. present day) 
gas and stellar surface density in the solar neighbourhood 
to their measured values (we assume T,gas ~ 14.5 Mo/pc^ 
(Oiling & Merrifield 1998), E.tars = 40 Mg/pc^ (Gould et al 
1996)). Figure 1^ illustrates several notable features. Firstly, 
different star formation prescriptions lead to very similar ex- 
ponential behaviours (over at least 5 scale lengths) for the 
stellar surface density profiles; models really only deviate 



from each other at early times and at large (unobservably 
faint) radii. Interestingly, the two component model pro- 
duces rather smooth profiles despite the abrupt change in 
the form of the star formation and viscosity laws at roughly 
three scalelengths. These particular profiles deviate some- 
what from exponential behaviour beyond 5-6 scale lengths, 
but the effect is subtle. Gas profiles are more sensitive to 
the different star formation laws, but generally exhibit ex- 
ponential declines in their outer regions with e-folding scales 
significantly larger than that of the stars. At late times, the 
scale length of the gas at large radii, where the gas is still 
the majority component, exceeds the stellar scale length by a 
factor of approximately 2. Secondly, we find the stellar com- 
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ponent exhibits self-similar growth out to ~ 5 scale lengths 
over the entire star forming period, implying an invariance in 
the stellar scale length with time, as first noted by Olivier et 
al (1991). This is further quantified in Figure ^ which shows 
that the scale length lies well within ~ 20% of its present- 
day value throughout the disc's history. Thus, while viscous 
discs are characterised by inside-out growth, which has led 
to the speculation that the resulting scale length should in- 
crease with time (eg. Zhang & Wyse 2000), we here show 
that this growth is simply manifest as an expansion in the ra- 
dial extent of the exponential stellar disc due to the viscous 
transport of gas from smaller radii. This effect is unlikely to 
be verifiable by studies of high redshift galaxies, however, 
since the strongest evolution occurs in regions of the disc so 
faint that they are inaccesible to direct observation. 

Figure ^ shows the star formation rate (dashed lines) 
normalised to the peak rate as a function of time at 2.7 ('so- 
lar neighbourhood'), 5 and 10 exponential scale lengths for 
the particular star formation law T.sfr oc Tigas/R; other 
star formation prescriptions exhibit similar behaviours. For 
a Freeman exponential disc (Es = 21.65 mag/D ), 5 scale 
lengths corresponds to an expected surface brightness of 
~ 1% of sky, or about the current limit for quantitative stud- 
ies of surface brightness in external galaxies. The top axis 
shows age of the disc in Gyr, where we have assumed the 
present-day stellar disc age of the Milky Way to be 11 Gyr as 
determined from the ages of the oldest stars in the Hipparcos 
sample (Binney, Dehnen & Bertelh 2000). The disc at 2.7 
and 5 scale lengths exhibits slowly declining star formation 
rates in the pre-assembled model, decreasing by factors of 
roughly 5 and 2 over the disc's lifetime respectively. On the 
other hand, the disc at 10 exponential scalelengths is char- 
acterised by a slowly increasing star formation rate (bottom 
panel Figure ^ reflecting the finite time required to trans- 
port gas outwards to these parts by viscous torques; the 
arrival of material in the extreme outermost disc is there- 
fore delayed long after the onset of star formation at smaller 
radii. This delay is not a major effect at 5 scale lengths 
however, even for the most compact plausible initial gas dis- 
tributions. 

Also shown in Figure ^ is the fraction of the final stellar 
mass in place as a function of time (solid lines) . These curves 
are contrasted with the analytic predictions of a simple tn 
situ model (dashed-dotted line) obeying the same star for- 
mation law (note that in the in situ model the total - gas 
-I- stars - radial profile is exponential by construction at all 
times) calculated as: 



dt 



R 



(3) 



where the adopted star formation efficiency, k, is chosen to 
be consistent with the viscous model. Integration yields 



'kt/R 



^stars i^tnow^ 



(4) 



where we have normalised the mass in stars present at a 
given time to that of the present day. While the mean age of 
the stellar population present at 5 scale lengths is somewhat 
smaller than at 2.7 scale lengths, the effect is rather small 
and the bulk of the stars at this radius have formed fairly 
early on. Specifically, while 50% of the final stellar surface 
density is in place after 0.3 of the star formation duration (or 
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Figure 3. Evolution of the stellar exponential scale length rela- 
tive to the present-day value as a function of time since the onset 
of star formation for the case of a pre-assembled gas disc. Results 
for different star formation laws are shown. At each epoch, the 
stellar scale length is measured by fitting the surface density pro- 
file over the region corresponding to ~ 0.3 — 0.9 the radius of the 
initial gas disc. 



~ 3.5 Gyr for our assumed disc age) at 2.7 scale lengths, this 
same amount is in place after ~ 0.4 of the star formation pe- 
riod (or ~ 5 Gyr) at 5 scale lengths. Interestingly, there does 
not appear to be any great difference between the growth 
rate of stellar mass at 5 scale lengths in the viscous and in 
situ models, implying that viscous fiows have little impact 
on the star formation history at these radii; observations of 
stellar age distributions alone in these parts could not dis- 
tinguish between viscous and in situ scenarios. Beyond the 
radius of the initial gas disc, the stellar population becomes 
increasingly young with radius. 

3.2 The Case of Extreme Infall 

We now contrast the models above with viscous models in 
which tin fall '~ t„om- Wc refer to this as the 'extreme in- 
fall' case, since the total mass infall rate onto the galaxy 
has been constant up to nearly the present day, and star 
formation in the disc started early on in the infall period 

{tstar TO^* 

nfaii)- Figure H shows the stellar and gaseous 
surface density profiles for different star formation prescrip- 
tions at various times since the onset of star formation. For 
all models, the final stellar and gaseous profiles are almost 
identical for the extreme infall viscous model and the pre- 
assembled viscous model (contrast the highest stellar profiles 
and lowest central gas profiles in Figures |2| and ^) . The tem- 
poral evolution of these profiles is very different however. 
While the disc is able to continuously adjust to the chang- 
ing profile of the infallen material so as to generate an ap- 
proximately exponential stellar profile at all times, the scale 
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Figure 4. The star formation history (dashed line) of the disc 
at radii corresponding to 2.7 (top), 5 (middle) and 10 (bottom) 
exponential scale lengths in the case of a pre-assembled gas disc. 
The star formation rate is normalised to the peak rate at each ra- 
dius. These models are for a star formation law T:sfr = ^gas/R- 
The top axis shows the age of the disc, assuming the present-day 
age of 11 Gyr (Binney et al 2000). Also shown is the fraction of 
the final stellar mass in place (solid line) as a function of time 
and the corresponding analytic prediction of an in situ model 
(dashed-dotted line) obeying the same star formation law. We do 
not show the in situ predictions for the disc at 10 scale lengths, 
since this lies well beyond the extent of the initial gas disc where 
gas would never reach in the absence of viscous flows. 



length of the exponential evolves systematically with time, 
simply reflecting the mean specific angular momentum of the 
material which has fallen in by that epoch. The stellar pro- 
flies are always exponential over at least 4-5 scale lengths. 
For the density and angular momentum distributions in the 
halo assumed here (p oc r~^,Q oc r~^), the total angular 
momentum scales with the square of the infallen mass, so 
that the mean speciflc angular momentum is simply propor- 
tional to the infallen mass. Our default infall prescription, 
which makes speciflc assumptions about the halo properties 



and the cosmological model (see Sec , implies a constant 
infall rate for t < tinfaii and hence the mean specific angular 
momentum of the infallen material rises linearly with time 
over this period. Figure ^ demonstrates that the stellar scale 
length indeed displays this behaviour. 

A fundamental assumption in our infall model is that 
the mean specific angular momentum of the accreted mate- 
rial systematically increases with time. We feel this is a valid 
assumption, regardless of the exact nature of the infall, and 
should be preserved. For example, if the dominant mode of 
infall is high velocity cloud accretion (Blitz et al 1999), one 
would still expect those clouds falling in at late times to have 
originated further from the Milky Way and to have suffered 
on average more tidal torquing than clouds which fell in at 
early times. While individual clouds will carry a range of 
angular momenta (and indeed in our models a given shell 
maps onto a range of radii), the mean trend should be for 
the average angular momentum to increase with time. It is 
worthwhile investigating how much of the scale length evo- 
lution is due to other aspects of our adopted infall model 
however. Figure |^ shows the growth of the stellar exponen- 
tial scale length for an isothermal halo model in which the 
total infall rate is allowed to decrease exponentially with 
time. In this case, the amount of recent infall (ie. t ~ tnom) 
is reduced relative to that in the constant infall case by a 
factor ~ 11 Gyr /Tin J, where 11 Gyr is the present-age of 
the disc and ri„/ the e-folding time of the infall rate. In the 
limit of Tinf » 11 Gyr, this infall prescription mimics the 
constant infall case whereas for Ti„/ << 11 Gyr it mimics 
the pre-assembled disc case, with essentially all the infall oc- 
curing early on; Figure ^ shows an intermediate case where 
Tinf =5 Gyr. We also examine an infall model in which the 
baryons in the halo are described by a Navarro- Frenk- White 
(NFW) profile in which p « (r/rs)"^(l + r/rs)~'^ with the 
parameter r^ defining the scale radius. Navarro, Frenk & 
White (1997) argue on the basis of high resolution N-body 
simulations that this form better characterizes the equilib- 
rium density profiles of dark matter halos than an isother- 
mal sphere. In terms of infall rates, it leads to a somewhat 
higher infall rate than the isothermal sphere at early times 
and somewhat lower at late times. Figure ^ indicates that 
the growth of the scale length in these alternate infall models 
is very similar to that in the constant infall case and thus 
that our results on the scale length evolution are robust. 
Contrasting Figures ^ and ^ with Figure ^ one concludes 
that if clear evidence were to emerge for an increase in stel- 
lar scale length of galaxies over cosmological time, it would, 
within the framework of our model, be an unambiguous in- 
dicator of the late infall of high angular momentum material 
and could not merely be attributed to the effects of viscous 
evolution. 

A consistent picture emerges from analysis of the star 
formation history at different locations in the disc. Figure 
^ shows the star formation rate (dashed line) normalised to 
the peak rate as a function of time for our default star for- 
mation law. Very few stars form initially at the solar neigh- 
bourhood as a direct result of the low gas density (only gas 
viscously transported from smaller radii is present). As the 
peak infall rate at this radius is reached (indicated by the 
arrow in the top panel of Figure the star formation rate 
rapidly increases and then levels off to an almost constant 
value until the present epoch. The fraction of the final stellar 
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Figure 5. Snapshots of the gaseous (left) and stcUar (right) surface density distributions at various fractional times (1/3, 2/3 and 1) of 
the star formation duration (which terminates at tnow) for the extreme infall model. In this model star formation starts shortly after 
the infall phase begins (tstar ~ jQtinfall- The radial coordinate is normalised to the stellar exponential scale length measured at the 
end of the model. The final solutions are those with the highest stellar profile at all radii and the lowest inner gas profile. Different star 
formation laws are shown in top, middle and bottom panels. 



mass in place as a function of time for singular isothermal 
sphere and NFW halo baryon distributions is also shown 
(solid and dashed-dotted lines respectively). As expected, 
the disc at 2.7 scale lengths is predominantly young in the 
extreme infall model, with only a small dependence on the 
specific form of the infall prescription. Specifically, 50% of 
the final stellar mass in place after 0.7 of the star forming 
duration, or ~ 8 Gyr if a disc age of 11 Gyr is assumed 
(note that in this case very few 11 Gyr stars are expected 
to exist outside the innermost regions of the disc, where the 
first infall is received). At 5 and 10 scale lengths, the stel- 
lar population is conspicuously young in the extreme infall 
model. The disc at 5 scale lengths receives material from 



direct infall of halo gas only at the very end of the infall 
phase, so that prior to this, gas only reaches this location as 
a result of outward viscous transport from an initially much 
more compact gas distribution. This is especially true of the 
disc at 10 scale lengths, which never receives direct infall in 
our model. Consequently, the dominant star-forming phase 
in these regions only gets under way at late times. We draw 
attention to the fact that this combined infall/ viscous model 
may produce an unusual chemical signature in the very outer 
regions, comprising an older metal rich population (formed 
from prc-enriched gas transported outwards from small radii 
at early times) and a young metal poor population contain- 
ing stars formed from recent unenriched infall. We plan to 



© 0000 RAS, MNRAS 000, 000-000 



The Evolution of Stellar Exponential Discs 9 




Figure 6. Evolution of the stellar exponential scale length rela- 
tive to the present-day value as a function of the star formation 
duration for the constant infall model. Results for different star 
formation laws are shown. At each epoch, the stellar scale length 
is measured by fitting the surface density profile over the region 
corresponding to 0.3—0.9 of the instantaneous infall radius. 

investigate detailed predictions for the chemical signatures 
of viscous evolution and infall in a future paper. 



4 DISCUSSION 

We have explored models for the evolution of disc galax- 
ies subject to star formation, viscous redistribution of an- 
gular momentum, and cosmologically-motivated gaseous in- 
fall. We have focused attention on the extreme cases of these 
models and the rather different predictions which arise for 
the temporal evolution of the exponential scale length, and 
the radial growth rate of galactic stellar discs. In this sec- 
tion, we discuss how these model predictions compare with 
relevant observational data for both local and distant disc 
galaxies. 

4.1 Constraints from the Local Universe 

Stellar populations in the Milky Way and other nearby spi- 
rals provide the 'fossil record' of star formation and galaxy 
assembly over cosmic lifetimes. The wealth of high quality 
data available for stars in our local neighbourhood of the 
Milky Way disc makes it a logical starting point for model 
comparisons. 

An important prediction of our models is the star for- 
mation history - which is reflected in the stellar age distribu- 
tion - at various radii. The pre-assembled and extreme infall 
models make very different predictions for the star formation 
history of the disc at 2.7 exponential scale lengths. In the 
former case, the star formation rate declines smoothly by a 



Figure 7. Evolution of the stellar exponential scale length rel- 
ative to the present-day value as a function of different infall 
prescriptions (ie. halo profiles). As in Figure joj infall has been 
occuring in these models since the onset of star formation. At 
each epoch, the stellar scale length is measured by fitting the sur- 
face density profile over the region corresponding to 0.3-0.9 of the 
instantaneous infall radius. 

factor of ~ 5 over the disc's lifetime with the consequence 
that a significant fraction of disc stars form early on. In the 
latter case, however, very few stars form initially with the 
rate then increasing steeply to a roughly constant value. As 
a consequence, most stars form late in this model. Quantita- 
tively, for an assumed disc age of 11 Gyr, more than 50% of 
the final stellar mass is in place at solar neighbourhood after 
~ 3.5 Gyr in the pre-assembled model whereas less than 5% 
is in place by same time in extreme infall model. 

It is well-established that stars with a considerable 
range in ages exist in our local neighbourhood, a non- 
negligible fraction of which appear to be 'old', ie. ages 
> 8 Gyr (eg. Edvardsson et al 1993, Binney et al 2000, 
Rocha-Pinto et al 2000). Rocha-Pinto et al (2000) have re- 
cently re-determined the star formation history of the local 
Galactic disk from a sample of 552 late-type dwarfs exhibit- 
ing chromospheric activity. Stellar ages are derived using 
a new metallicity-dependent chromospheric activity-age re- 
lation (see their Table 3). Assuming this sample provides 
an unbiased representation of the total thin disc popula- 
tion, we use their age distribution to infer the cumulative 
mass in stars present in the solar neighbourhood as a func- 
tion of time as shown in Figure ^. This curve shows a good 
agreement with the predictions of our extreme infall viscous 
model. While the Rocha-Pinto data suggest a slightly higher 
star formation rate at early times (hence a small excess of 
very old stars), the overall similarity of the two star forma- 
tion histories at later times is quite striking. We reiterate 
that our star formation history is largely determined by the 
infall history at this radius, and not by viscous flows. Al- 
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Figure 8. The star formation history (dashed Une) of the disc 
at radii corresponding to 2.7 (top), 5 (middle) and 10 (bottom) 
exponential scale lengths in the extreme infall case. The star for- 
mation rate is normalised to the peak rate at each radius. These 
models are for a star formation law Y^sFR = ^gas/R- The top 
axis shows the age of the disc, assuming the present-day age of 
11 Gyr (Binney et al 2000). Also shown is the fraction of the final 
stellar mass in place for the particular cases of baryonic infall from 
halos with isothermal and NFW forms (solid and dashed-dotted 
lines respectively) as a function of time. The time corresponding 
to the peak infall rate at the solar neighbourhood is marked by 
the arrow in the top panel. 



though the possibihty remains that some observational bi- 
ases exist within the Rocha-Pinto et al sample - for example, 
the onset and duration of stellar chromospheric activity is 
poorly understood (eg. Balnius & Vaughan 1985) - it seems 
unlikely that these biases could alter the derived star for- 
mation history enough to make it more closely resemble the 
predictions of our pre-assembled model. 

Knowledge of the age distributions of stars at large radii 
in discs (for example, at 5 scale lengths) would have even 
greater leverage for model comparisons, since it is precisely 
these parts where the predictions of the extreme infall and 
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Figure 9. The cumulative mass in stars as a function of time in 
the solar neighbourhood as derived from the chromospheric age 
distribution of Rocha-Pinto et al (2000) (dashed line). Also shown 
are the predictions for our extreme infall model (dashed-dotted 
line) and the pre-assembled disc model (solid line). 



pre-assembled models deviate the most. Unfortunately, es- 
sentially no information is currently available on the range 
of stellar ages present in these parts. Due to our location 
within the Milky Way disc, this issue may be more easily 
addressed in the future by deep studies of the resolved stel- 
lar populations in nearby external galaxies than in our own. 

The effects of angular momentum transport alone mean 
that viscous models naturally predict inside-out growth for 
galactic discs whether or not infall is also present. Associat- 
ing the 'mean age' of the stellar population with the age by 
which 50% of the final stellar mass is in place, our models 
suggest an age gradient of ~ —0.7 Gyr over the radial raiige 
spanning 2-5 scale lengths in discs (see Figures ^ and ||). 
Inside-out formation has long been suspected in disc galax- 
ies on the basis of radial colour and metallicity gradients (eg. 
Tinsley & Larson 1978). For example, de Jong (1996b) finds 
that radial colour gradients (in the sense of bluer colours at 
larger radii) extending over several scale lengths are common 
in spirals and are largely consistent with trends in stellar age 
alone. Using stellar population synthesis models. Bell and de 
Jong (2000) have explicitly quantified the radial age varia- 
tion implied by these colour gradients and find the typical 
large spiral to exhibit an age gradient of ^ —0.8 Gyr per 
K-band scalelength, which also agrees with the predictions 
of our models. 

Finally, we note that there is still significant debate as 
to the importance of continued gaseous infall for disc galaxy 
evolution, with few existing direct observational constraints. 
Galactic high velocity clouds may be a signature of gas ac- 
cretion onto the Milky Way, but the lack of information 
concerning the distances of these clouds hinders a clear un- 
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derstanding of their origin. Blitz et al (1999) have carried 
out a simple simulation of the formation and evolution of 
the Local Group in which 10® particles (each representing 
a 10® Mq gas cloud) are subject to the gravitational forces 
of M31, the Milky Way and the tidal field of the other Lo- 
cal Group members. Any cloud falling towards the Milky 
Way or M31 which passes within 100 comoving kpc of their 
centres is assumed to be accreted by that galaxy. They find 
present-day neutral gas accretion rates of ~ 1 M0 yr~^ for 
both galaxies, with the rates being somewhat higher in the 
distant past. The overall evolution is consistent with an ex- 
ponential decline with an e-folding time of about 5 x 10^ yr. 
On the other hand, our extreme infall model, when nor- 
malised to the Milky Way, requires a constant infall rate of 
4.5 Mq yr^^. Our infall models therefore require gas accre- 
tion rates that are slightly higher but of the same order of 
magnitude as those found in Blitz's simulation. 

4.2 Constraints from Distant Galaxies 

An exciting recent development has been the availability 
of quantitative structural parameters for moderate redshift 
disc galaxies which can be used to directly map the size 
evolution of galaxies with cosmic epoch (eg. Roche et al 
1998, Lilly et al 1998, Simard et al 1999, Giallongo et al 
2000). Based on a sample of discs drawn from the CFRS 
and LDSS redshift surveys, Lilly et al (1998) find the size 
function of disc scale lengths in disc-dominated galaxies (ie. 
bulge-to-total ratios < 0.5) to be roughly constant to z~ 1, 
at least for large discs with R_d > kpc. In a similar 

study, Simard et al (1999) analyse the size-magnitude rela- 
tion defined by disc-dominated field galaxies in the DEEP 
survey. Carefully considering possible selection effects, these 
authors also conclude that there is no significant evidence 
for any evolution in the size-magnitude distribution of disc 
galaxies over the redshift range 0.1 < z < 1.1 and that a sig- 
nificant number of discs lie close to the canonical Freeman 
relation at all redshifts probed. 

Taken together, these studies suggest only a very mod- 
est change in stellar exponential scale length with epoch out 
to z ~ 1 (however see Mao et al (1998) for an alternative 
interpretation of the same dataset). For an Einstein de Sit- 
ter Universe, the interval between a redshift of unity and 
the present epoch can be associated with the period from 
~ 0.3 — 1.0 times the star formation interval in our models, 
assuming that star formation in the disc begins early on. 
Inspection of Figures ^ and ^ indicates that the disc scale 
lengths are expected to increase by about a factor of 3 over 
this redshift range if gaseous infall is important. Thus, in 
contrast to the Milky Way constraints discussed in the pre- 
vious section, it would appear that observations of moderate 
redshift galaxies lend support for a scenario in which con- 
tinued infall - as we have modelled it - is not important for 
building galactic discs. We note that the global star forma- 
tion rates of the discs which form in our pre-assembled and 
extreme infall models decrease and remain roughly constant 
with time respectively, neither of which are inconsistent with 
current observations of large disc galaxies. 

It should be kept in mind, though, that quantitative 
characterization of distant galaxy profiles is still plagued 
by a number of difficulties. For the highest redshift galax- 
ies in the aforementioned samples, exponential profiles are 



generally fit over 1-2 disc scale lengths and hence are par- 
ticularly susceptible to errors in the bulge and disc separa- 
tion. Furthermore, studies to date have measured structural 
parameters in a single passband over the redshift range of 
interest. If strong colour gradients are present in these sys- 
tems, this could conspire with the fc-correction to mimic a 
constant scale length in a fixed passband over a range in red- 
shift even though the rest-frame scale lengths are evolving. 
Colour gradients significantly in excess of those measured in 
local galaxies would be required for this to be an important 
effect however. 



5 CONCLUSIONS 

We have explored models for the evolution of disc galaxies 
that invoke simultaneous star formation and viscous redistri- 
bution of angular momentum, as well as a cosmologically- 
motivated prescription for gaseous infall. The existence of 
galaxies which exhibit exponential disc profiles over signifi- 
cantly more than ^ 5 scale lengths (eg. Weiner et al 2000, 
Barton & Thompson 1997) lend strong support for viscous 
models as viscosity may be the only viable mechanism for 
transporting material to these extreme radii, while maintain- 
ing a smooth exponential profile in the stellar disc. Following 
Lin & Pringle (1987a), we fix the relationship between the 
star formation and viscosity prescriptions such that the char- 
acteristic timescale for material to move a fractional radius 
of order unity is comparable with the timescale on which it 
is converted into stars (i.e. ^ t^). While the inclusion of 
viscous fiows is essential for ensuring that the stellar disc is 
always exponential, we find that such flows play essentially 
no role in determining the evolution of the disc scale length, 
or the star formation history over the inner ~ 5 scale lengths 
(corresponding to the maximum infall radius in our model) . 
Indeed, the dominant process governing the size evolution 
and star formation history of the luminous stellar disc is 
that of gaseous infall. 

We summarize our main results as follows: 

1. In models in which the main infall phase precedes the on- 
set of star formation (ie. a pre-assembled gas disc), we find 
the exponential scale length to be rather invariant with time 
over the entire history of the disc. This result holds true for 
a variety of popular star formation prescriptions. The value 
of the stellar disc scale length is primarily determined by 
mean specific angular momentum of the initial gas disc. 

2. In the pre-assembled model, the star formation rate 
smoothly declines with time within the inner 5 scale lengths 
(approximately the maximum infall radius, Vinit)- As a re- 
sult, there is a considerable 'old' (ie. > 8 Gyr) population 
present at these radii, with half of the final stellar mass in 
place more than 6 Gyr ago. The onset of star formation is 
delayed at radii greater than Vinu as the outer disc is formed 
entirely by material viscously transported from smaller radii; 
the stellar populations in these parts are thus increasingly 
young with radius. 

3. If infall continues during the main star-forming phase 
of the galaxy's evolution, the stellar disc remains approx- 
imately exponential at all times, but now the scale length 
increases towards the present epoch, in simple proportion to 
the average specific angular momentum of the material fallen 
in to the galactic disc to date. For example, for an isothermal 
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halo density profile and constant spin parameter, the scale 
length is directly proportional to the infallen mass which is 
proportional to time. Regardless of whether viscosity is at 
work or not, late infall of high angular momentum material 
is the only way to increase the stellar scale length of discs 
within the framework of our models. 4. In contrast to the 
case of a pre-assembled disc, the disc stellar populations in 
the infall models are predominantly young, with more than 
half of the stars even at the location of solar neighbourhood 
predicted to have ages < 5 Gyr. The disc at 5 scale lengths 
receives direct infall of halo gas only at the very end of the 
accretion phase, so that prior to this, gas only reaches this 
location as a result of outward viscous transport. As in the 
pre-assembled case, all material at radii beyond the maxi- 
mum infall radius has been viscously transported there from 
smaller radii. 

5. We compare our model predictions with recent observa- 
tions of disc galaxies sizes and stellar populations. A par- 
ticularly good agreement is found between the solar neigh- 
bourhood star formation history predicted by our extreme 
infall model and the recent observational determination of 
this quantity by Rocha-Pinto et al (2000). On the other 
hand, the relative invariance of disc scale length with cosmic 
epoch, as derived from studies of moderate redshift galaxies, 
appears to lend support for our pre-assembled model. Fu- 
ture observations that will be important for clarifying this 
situation include detailed studies of resolved stellar popu- 
lations at large galactocentric radii (^5 scale lengths) in 
the Milky Way and other nearby spirals in order to con- 
strain the age and star formation history of the outer disc, 
and larger, deeper surveys of the structural parameters of 
moderate redshift disc systems in various passbands. 
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